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Abstract 

We show the dynamical stability of a six-dimensional braneworld solu¬ 
tion with warped flux compactification recently found by the authors. 
We consider linear perturbations around this background spacetime, as¬ 
suming the axisymmetry in the extra dimensions. The perturbations 
are expanded by scalar-, vector- and tensor-type harmonics of the four¬ 
dimensional Minkoswki spacetime and we analyze each type separately. It 
is found that there is no unstable mode in each sector and that there are 
zero modes only in the tensor sector, corresponding to the four-dimensional 
gravitons. We also obtain the first few Kaluza-Klcin modes in each sector. 


1 Introduction 


It is being established that there are two major accelerating, or quasi de Sitter, phases 
in the history of our universe. One is inflation in the early universe and the other is 
the current phase of accelerating expansion. This picture is supported by observations 
such as cosmic microwave background and supernovae dim and thus we are almost 
sure that there are and were such accelerating phases. However, we essentially do not 
know what causes those phases, i.e. we do not yet know what inflaton and dark 
energy are from the viewpoint of fundamental physics. 

In string theory it had been thought rather difficult to construct a four-dimensional 
de Sitter universe with stabilized moduli until the construction by Kachru, Kallosh, 
Linde and Trivedi (KKLT) |(3j. (See 0, (i for followup proposals.) They evaded pre¬ 
viously known no-go theorems by putting anti-D-branes at the bottom of a warped 
throat after stabilizing all moduli. Since the shape of the warped, compact extra 
dimensions are stabilized by fluxes, this set-up is often called warped flux com- 
pactihcation. This set-up provides a number of possible applications to cosmol¬ 
ogy 0 Em imiiii]. 

In our previous paper El we pointed out that brane gravity in the warped flux 
compactihcation is somehow similar to that in the first Randall-Sundrum (RSI) sce¬ 
nario E with radion stabilization [TT] . (See refs. JT3I TS. (T71 -151 for brane gravity 
in the RSI scenario with radion stabilization.) The evolution of matter on the brane 
changes the bulk geometry not only near the brane but possibly everywhere in the 
whole extra dimensions. Provided that all moduli are stabilized, the bulk geometry 
should quickly settle to a configuration which is determined by the boundary con¬ 
dition, i.e. the brane source(s), values of conserved quantities and the regularity of 
the other region of the extra dimensions. As a consequence of the change of the 
bulk geometry, the induced geometry on the brane responds to the evolution of the 
matter source on the brane. The four-dimensional Einstein theory is recovered as a 
rather indirect and subtle relation between the matter source on the brane and the 
response of the induced geometry. To support this picture, we considered a simplified 
situation in which we can see the recovery of the four-dimensional Einstein theory in 
the warped flux compactihcation. In particular, we found an exact solution repre¬ 
senting a six-dimensional brane world with warped flux compactihcation, including a 
warped geometry, compactihcation, a magnetic hux, and one or two 3-brane(s). (See 
refs. [HE EH EH 1221231 El for other models of six-dimensional brane world.) 

The purpose of this paper is to show the stability of the six-dimensional exact 
solution found in the previous paper [T2] , For simplicity we set the four-dimensional 
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cosmological constant to zero and assume axisymmetry in the extra dimensions. We 
expand linear perturbations by scalar-, vector- and tensor-type harmonics of the four¬ 
dimensional Minkowski spacetime and analyze each type separately. Linear pertur¬ 
bations in six-dimensional models are also considered in [25112 6| . We shall show that 
there is no unstable mode in each sector and that there are zero modes only in the 
tensor sector, corresponding to the four-dimensional gravitons. We also obtain the 
first few Kaluza-Klein modes in each sector. 

The rest of this paper is organized as follows. In section [21 we briefly review our 
six-dimensional brane world model. We then consider linear perturbations with ax¬ 
isymmetry to show the stability of our exact solution in section [21 We numerically 
solve the perturbed Einstein equations and Maxwell equations for each type of per¬ 
turbations, and then show that there is no unstable mode. In sectional we summarize 
the main results and discuss them. 


2 6D Warped Flux Compactification 


In this section, we briefly review our six-dimensional braneworld model, which cap¬ 
tures essential features of the warped flux compactification, that is, warped geome¬ 
try, magnetic flux of an antisymmetric held along the extra dimensions, and branes. 
The model is simple enough to make it possible for us to analyze gravity on the 
brane from higher dimensional point of view. Our start point is the six-dimensional 
Einstein-Maxwell action, 



( 2 . 1 ) 


where M 6 is the six-dimensional reduced Planck mass, A 6 is the bulk cosmological 
constant, and Fmn = 9m An — 9nAm is the held strength associated with the 1/(1) 
gauge held Am- 

The bulk solution considered as the background in this study is 


/ / / ■ - 

ds\ = r 2 r)^dx fJ ’dx 1 ' + —— + f(r)dcj) 2 , 

f ( r ) 

A M dx M = A(r)d(f>, 


( 2 . 2 ) 

(2.3) 


where 
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This solution corresponds to the A 4 ± —» 0 limit of the more general solution found 
in |Q, where A 4 ± is the four-dimensional cosmological constant on the branes. (See 
Appendix A.2. of |Q for this limit.) As explained in P3, this solution is related by 
a double Wick rotation to a topological black hole with the horizon topology R 1 . 

Next, we consider embedding of one or two 3-brane sources. As is well known, an 
object with codimension 2 induces a deficit angle around it as 


5 


<7± 

ur 


( 2 . 6 ) 


where 5± is the deficit angle due to the tension o± of the brane. We should note 
that this formula is valid under the axisymmetry if radial stress is much smaller than 
energy density m- 

Hereafter, we assume that the function /(r) given by Eq. (12. 4 j) have two positive 
roots r = r± (0 < r_ < r + ) for f(r) = 0 and is positive between them (r_ < r < r+). 
This requires that the six-dimensional cosmological constant A 6 be positive, and thus 
we assume A 6 > 0 throughout the paper. Since / vanishes at r±, r — r± defines 
surfaces of codimension 2. Thus we can put the 3-branes at r±. The periods of the 
angular coordinate A0 calculated at r± must coincide, thus we get 


A(j) 


2n — S. 
«+ 


where <5± are the deficit angles, 


27T — <5_ 


(2.7) 


K ± = Tj/Vi), 

and a prime denotes derivative with respect to r. This is rewritten as 

27T — <5 + n + 

2tt — K ._ ’ 


( 2 . 8 ) 


(2.9) 


which can be regarded as a boundary condition since the l.li.s is specified by the brane 
sources and the r.h.s. can be written in terms of the bulk parameters fib and b. 

The background geometry can be specified by the three parameters, the tensions 
cr± of the brane at r± and the magnetic flux $ in the bulk. The magnetic flux $ is 
given by the equation 


$ 

-v/(2jr - <5 + )(2>r - S_) 


_ i _ (l 



( 2 . 10 ) 


where L = ^10/Ag and the r.h.s is written only by fib and b. If we specify a± and 
<f), the left hand sides of Eq. m and Eq. are determined. On the other hand, 
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the right hand sides of these equations are written by pb and b. Thus we can solve 
(EH) and (12.101) w.r.t. pb and b. We can also determine A0 by Eq. (12. 7[) . 

Before going into details, we make some general remarks on the background solu¬ 
tion. Our exact braneworld solution captures some essential features of the warped 
flux compactification, including warped geometry, compactification, moduli stabiliza¬ 
tion, flux and branes. However, since the higher-dimensional cosmological constant Ag 
is assumed to be positive, we have to confess that this model of warped flux compacti¬ 
fication does not completely mimic the KKLT construction, where the 10-dimensional 
cosmological constant is not positive but zero. Nonetheless, our model is useful in the 
sense that it provides a testing ground on which brane gravity with warped flux com¬ 
pactification can be analyzed from higher dimensional viewpoints. Note that many 
cosmological considerations in the KKLT setup are based on the implicit assumption 
that 4-dimensional Einstein gravity should be recovered at low energy and that it is 
very important to see whether this assumption is viable or not from higher dimen¬ 
sional viewpoint. As discussed in D2. the recovery of the 4-dimensional Einstein 
gravity in warped flux compactification is not as simple as would be expected from 
the 4-dimensional effective theory. The 4-dimensional gravity is recovered in a rather 
subtle way as a consequence of bulk dynamics, and the exact solutions in our model 
were useful to see how this picture works explicitly. 

Finally, we explain the limit a = r_/r + —> 1 of the bulk geometry. In this limit, 
the coordinate distance r + — r_ between the two branes vanishes, and thus the bulk 
geometry appears to collapse. However, the proper distance between r = r_ and r + 
does not vanish and the geometry of extra dimensions remains regular. This can be 
explicitly shown by a coordinate transformation: 


2 r — (r + + r_) 

r = --- 

r + — r_ 

ip = A 6 (r+ — r-)(j> 


(-1 < r < 1), 


( 2 . 11 ) 

( 2 . 12 ) 


Even if we take a —> 1, the domain of r remains finite while that of r vanishes. With 
this new coordinate system, the metric of the extra dimensions becomes 



(2.13) 


where 


f(r 


f(r) 
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and we defined 


P-r + P+ 

2 


2 . 


(2.14) 


P± = or 1/2 ±a 1/2 , (2.15) 

2 n 

7n = ^y-". ( 2 - 16 ) 

i=0 

The function / includes only the background parameter a , and is defined so as to 
remain finite in the limit of a —► 1. Indeed, / in this limit is as simple as 

/=y(l-r 2 )- (2-17) 

Thus, the metric $nnp becomes that of a round sphere of radius 1/a/2A 6 in this 
limit. 

As we mentioned, the coordinate f runs over the finite interval [—1,1]. On the 
other hand, the period of the coordinate <p appears to collapse since the coefficient 
(r + — r_) in the definition (EH2D vanishes in the a —> 1 limit. Actually, this is not 
the case. The “surface gravity” k± defined in Eq. (ESI) is written in terms of / as 

K± = T(r+-r_)dJ(± 1). (2.18) 

Therefore, the new coordinate ip has the period 

Av 3 = A 6 (r + - r_) A(j) = (2vr - <J±), (2.19) 

which is indeed finite, and becomes Atp = 2ir — <5+ = 2(/r — <5_ in the a —> 1 limit. 
Thus, the geometry of extra dimensions is nothing but a round sphere with a deficit 
angle 5 + = <5_, i.e. a football-shaped extra-dimensions considered in EDI EH- 

As we shall see in the next section, when the new coordinates (f, p >) are used and 
the KK mass is properly rescaled, the system of the perturbation equations depends 
on the background parameters through just one parameter a. Thus, we calculate the 
mass spectra of KK modes in each type of perturbations while changing a from 1 to 
0. The perturbation equations written in terms of (f, ip) remain regular in the a —» 1 
limit. 


3 Dynamical stability 

In this section, we consider linear perturbations around the background spacetime 
described in the previous section in order to show the dynamical stability of this 
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spacetime. For simplicity we assume that the perturbations are axisymmetric in the 
extra dimensions. Since the background geometry has the 4-dimensional Poincare 
symmetry, it is convenient to expand perturbations by scalar-, vector- and tensor- 
type harmonics of the 4-dimensional Minkowski spacetime. The perturbations are 
labeled by its type and values of mass squared m 2 = —ri lw k l ,k l ,, where k M is the 4D 
projected wave number of each harmonics. 

To study the stability, we take advantages of the fact that for any perturbation 
type, the perturbation equations are reduced to eigenvalue problems with eigenvalue 
m 2 . We regard the background spacetime to be dynamically stable if the spectrum 
of m 2 is real and non-negative. In each sector, our strategy to tackle this problem 
consists of the following four steps: 

(i) We show the reality of m 2 . 

(ii) We rewrite the systems of the perturbation equations into a form which depends 
on background parameters through just one parameter a. The parameter a runs 
over the finite interval [0,1]. 

(iii) In the a —> 1 limit we analytically solve the perturbation equations and show 
that the spectrum of m 2 is non-negative. 

(iv) We numerically evaluate how each eigenvalue m 2 changes as a changes from 1 
to 0, and show that the spectrum remains non-negative throughout. 

The steps (i) and (ii) make it possible for us to analyze the stability in the a-m 2 
plane, i.e. a 2-dimensional space spanned by the paramter a and the eigenvalue m 2 . 
In particular, each eigenvalue generates a curve in the a-m 2 plane as a runs over the 
interval [0,1], and what we have to show is that all such curves are in the stable region 
m 2 > 0. Thus, provided that each eigenvalue is a continuous function of a , the step 
(iii) implies that, if there exists a value of a (0 < a < 1) for which the background 
is unstable, the curve of the lowest eigenvalue must cross the a axis at least once. 
Finally, the step (iv) proves the stability of the system. In this way, we shall show that 
our six-dimensional brane world model is dynamically stable against tensor-, vector- 
and scalar-type linear perturbations. While we can easily show the stability against 
tensor-type perturbations even without solving the perturbed Einstein equation, we 
need to perform numerical calculations in the step (iv) for scalar- and vector-type 
perturbations. For completeness, we shall perform numerical calculations for tensor 
perturbations as well and show first few KK modes. 

In the analysis of each sector, we shall require regularity of physically relevant, 
geometrical quantities such as the Ricci scalar of the induced metric on the brane, 
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the tetrad components of the bulk Weyl tensor evaluated on the brane, etc. This 
is because we are adopting the thin brane approximation, i.e. (ESI), and all we can 
and should trust is what is obtained within the validity of this approximation. If e.g. 
the Ricci scalar of the induced metric were singular then our approximation would 
be invalidated. It is of course possible to regularize the singularity by introducing 
a finite thickness of the brane. However, in this case the natural cutoff of the low 
energy effective theory is the inverse of the thickness, and in general the regularized 
“would-be singularity” is not expected to be below the cutoff scale. This simply means 
that we need a UV completion, e.g. the microphysical description of the brane, to 
describe the physics of the regularized “would-be singularity”. Therefore, in general 
we have two options: (a) to specify a fundamental theory such as string theory as 
a UV completion and go on; or (b) to concentrate on modes which are within the 
validity of the effective theory. In this paper we adopt the latter attitude, assuming 
the existence of a good UV completion but never using its properties. This is the 
reason why we adopt the thin brane approximation and require the regularities. 

Let us now start the analysis. To begin with, we expand the perturbed metric by 
harmonics of the 4-dimensional Minkowski spacetime as 

SgMNdx M dx A = h rr Y dr 2 + 2 h r ^Ydrd(j) + h^Y dcp 2 

+ 2 [h{T) r V{T)n + h^L) r V(^L)i?jd:rdx^ + 2 + ^(L)</>U(L)/^ d(f)dx^ 

+ (h(T)T(T)fiv + h(LT)T(LT)nv + h(LL)T(LL)in> + d(Y)T(Y)fj,uj dx^dx”, (3.1) 

where Y, V(t,l) and T(t,lt,ll,y) are scalar, vector and tensor harmonics, respectively. 

See Appendix IA.1I for definitions of the harmonics. The coefficients h rr , h r $, h 
h{T,L)ri d(T,L)<t> and h(T,LT,LL,Y) are supposed to depend only on r. In the same manner, 
the perturbations of the U( 1) gauge field can be expanded as 

5A M dx M = a r Ydr + a^Ydcj) + (a( T )V( T ) M + a(L) (3.2) 

Here the coefficient a r , a$ and a(r,L) are also supposed to depend only on r due to 
axisymmetry. 

The Einstein equations and the Maxwell equations are decomposed into three 
groups, each of which contains only variables belonging to one of the following three 
sets of variables {/i(T)}, {h(r) r , a (T)} and { h rr , h r </,, h^, h^y, h( LLY ), a r , a^, ci(L)}- 

Variables belonging to each set are called tensor type, vector type and scalar type, 
respectively. It should be noted that these variables include degrees of freedom of 
gauge transformation, which are explicitly given in Appendix IA.2I 
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3.1 Tensor-type perturbation 


Our first task is to show that m 2 > 0 for any non-vanishing tensor perturbations. We 
first derive the evolution equation for tensor perturbations in subsection Id. 1.1 1 From 
this equation and regularity of h and h', we can show m 2 > 0 without solving the 
evolution equation. We also show in subsection Id. 1. II that the system of the evolution 
equation and the boundary condition can be rewritten in the form which depends on 
the background parameters through just one parameter a = r_/r + . Then, we seek 
the analytic solution for a = 1 in subsection Id. 1.21 Using this result, we numerically 
solve the perturbation equations by relaxation method in subsection Id.l.dl Detailed 
explanation of relaxation method is given in Appendix lA.dl 


3.1.1 Basic equations 


For tensor perturbations, 

dsl 

A M dx M 


r 2 {r]^ + hYs^)dx ti dx v + — + fd<j) 2 , 
Adq i, 


( 3 . 3 ) 


where the perturbation is specified by the function h of r and the harmonics Y = 
exp(f/c M aU). The symmetric polarization tensor satisfies = k^s^ v = 0 for 
kJ'kn ^ 0, or = k^s^ = = 0 for k^k^ = 0, where is a constant timelike 

vector. There is no relevant equation coming from the Maxwell equation, and the 
perturbed Einstein equation becomes 

( r 4 fh + m 2 h = 0, (3.4) 

where a prime denotes derivative with respect to r and m 2 = —r] fJ ’ u k^k u . The regu¬ 
larity of the four-dimensional Ricci scalar of the induced metric on the brane requires 
that h should be regular. With the above equation for h, the regularity of the tetrad 
components of the six-dimensional Weyl tensor on the brane requires that b! should 
also be regular. Indeed, provided that h is regular, C r ^ ru /r is regular at r = r± if and 
only if h' is regular. 

With the regularity of h and h' at r = r±, it is easy to show that m 2 > 0 for any 
non-vanishing solutions: 

m 2 f dr r 2 h 2 = — f dr h (■r 4 //^ / ) , = / dr ( r 2 h') 2 f > 0, (3.5) 
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where we have done integration by parts and used the fact that f(r±) = 0. The 
equality holds if and only if h' = 0 in the region r_ < r < r + . Thus, there is no 
instability in the tensor sector and the zero mode is h = const. 

The above differential equation (ED include the two background parameters 
(Hb,b), or equivalently (r + ,r_). From now, we show this equation can be rewritten 
into a form which includes only one parameter a = r_/r + . First of all, we perform a 
coordinate transformation Eq. (IXTTD . In terms of r, the Eq. (TTTll becomes 

lo _ 1 q V 2 d r ((P-r + P+f fdfh) + rh 2 h = 0, (3.6) 

(P-r + P+) 


where we dehned 


o 

9 m 

m = -. 

r + r_ 


(3.7) 


The function / and thus the Eq. (El include only the parameter a. The boundary 
conditions for h are obtained by assuming that h can be expanded in the Taylor series 
at r = r±: 


3fh T 


m 


{p-f+p+y^f 


-.h 


= 0. 


(3.8) 


+±i 


Thus, it is sufficient that we solve the perturbation equation with a variety of a. 


3.1.2 Analytic solution for a — l 

Here, we give the solution of h for a = 1, which is used when we numerically solve 
the perturbation for general a in the next subsection. Substituting / into Eq. ESD 
and taking a —> 1, we get 


df ((l — r 2 ) dfh) + /i 2 h = 0, 


(3.9) 


where 
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This can be solved analytically as 


h = CP U (f) + DQ U (f), 


(3.10) 


(3.11) 


where C and D are normalization constants. P v and Q v are the Legendre functions 
of the first and second kind, respectively, v is related to p 2 as p 2 = v(y + 1). The 


9 








regularity of h at f = 1 requires D = 0, and that at f = — 1 does P u not to diverge 
there. This is realized if only if v is non-negative integer. The mass spectrum is then 
determined as 

fj, 2 = i/(i/+ 1) (i/ = 0,1,2, (3.12) 

The zero mode (a mode with m 2 = 0) is h = const, and all other modes (i.e. Kaluza- 
Klcin modes) have positive m 2 . Tims, the background spacetime for a = 1 is dynam¬ 
ically stable against tensor perturbations. 

3.1.3 Numerical solution of KK modes 

Here we obtain the first few KK modes of tensor-type perturbation by numerically 
solving the perturbed Einstein equation (EH) with the boundary condition (EHD- For 
this purpose, the relaxation method is useful. Detailed explanation of the relaxation 
method is given in Appendix I A. 31 Here we give an outline of this method. We first 
rewrite the second order differential equation (EH) to a system of two first order 
differential equations by defining d f h as well as h as a dependent variable. Then, 
the differential equations are replaced by finite-difference equations on a mesh of M 
points that covers the range of the integration. We start with an arbitrary trial 
solution which does not necessarily satisfy the desired finite-difference equations, nor 
the required boundary conditions. The successive iteration, now called relaxation, 
will adjust all the values on the mesh so as to realize a closer agreement with finite- 
difference equations and, simultaneously, with the boundary conditions. Good initial 
guesses are the key of efficiency in the relaxation method. Here we have to solve the 
problem many times, each time with a slightly different value of a. In this case, the 
previous solution will be a good initial guess when a is changed, and it will make 
relaxation work well. As shown in the previous subsection, the perturbation equations 
can be analytically solved for a = 1. Thus, we solve the problem while changing a 
from 1 to 0. 

We show the first four KK mode solutions of h for a given value of a in figure Cl 
When we plot the solutions, normalization is determined by using the generalized 
Klein-Gordon norm. See Appendix I A. 41 for its derivation. For tensor perturbations, 
it is defined by 

($, T) kg = i ~~~~~~ I d 3 x J drr 2 rf l,J 'rr' - %, v d f ^c'A • (3-13) 

Substituting = f d^k h^s^Y and = f d A k' , we get 

^ )kg = ( k o + K) ( k - k 0 J drr 2 h ni (r)h n2 (r), (3.14) 
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where subscript rik means h nk is the solution of tensor perturbations with the eigen¬ 
value , and we normalized the symmetric polarization tensor as s fjv s ,iu = 1. In 
terms of the coordinate (f, tp), this becomes 


(*,*)«? = (k 0 + K)6 3 ( k-\d)Atp 


MqT + t_ f ((3_r + (3 + ) 


16A 6 

r4„ 


dr- 


IT 

= (fc 0 + fc')5 3 (k-k')Av9-^f^(/r ni ,h n2 ). 

16A 6 


We normalize the solution by 


(h-n 1 1 h n , 

(U) 


= s. 


711712 J 


-h ni (r)h n2 (r) 

(3.15) 

(3.16) 


so that the zero mode solution is normalized as h(r ) = 1. We can easily prove the 
orthogonality between modes with different m 2 by using the equation of motion for 
h, (13.6p . 

Finally, we show the spectrum of the mass squared of first four KK modes as a 
function of a. The mass squared that has physical meaning is m 2 ± = — r^rj^kukj,, 
which are the ones observed on the brane at r = r±. They are related to rh 2 as 
m± = a ±lr rh 2 . We plot the spectrum of m 2 + in £gure|21 As is easily seen, m 2 + is non¬ 
negative for the entire range of a. Therefore, the background spacetime we consider 
is dynamically stable in the tensor-type sector. 

We also notice that m 2 + remains finite in a —> 0 limit, which implies m 2 + oc a 0 , 
whereas m 2 _ cx a~ 2 for a —> 0. This result is consistent with our previous study Q. 
where we analyzed how the Hubble expansion rate H± on each brane changes when 
the tension of the brane changes. In that paper we also considered higher-order 
corrections of the effective Friedmann equation with respect to H±. The result is 
that the higher-order corrections appear when H± get larger than critical values H*±. 
For a —» 0, we found that H*± behave as 


H 2 : oc a 0 , 


H 2 oca" 2 . 


(3.17) 


Since the higher-order corrections are caused by the KK modes, this energy scale 
corresponds to their mass measured on each brane. Thus, the behavior of rri± we 
obtained in the a —> 0 is consistent with our previous result in 12). 


3.2 Vector-type perturbation 

Next we show that m 2 > 0 for any non-vanishing vector perturbations satisfying 
relevant boundary conditions. We first derive the perturbed Einstein equations in 


11 








Figure 2: The spectrum of m\ for tensor perturbations as a function of a. 

subsection 13.2.11 There are two physical degrees of freedom for vector perturbations. 
Their differential equations are easily written into a form which is manifestly hermite, 
therefore m 2 is real. They are also written into a form including only the background 
parameter a. We then derive the boundary conditions from those equations, by 
assuming that the two perturbation variables can be Taylor expanded with respect 
to f ± 1. In subsection we summarize the analytic solution for a = 1. Using 

this result, we numerically solve the perturbation equations by relaxation method in 
subsection 13,2.31 

3.2.1 Basic equations 

For vector perturbations, 

dv 2 

ds 2 6 = r 2 r}^ u dx^dx v + 2 {h^ydr + h^d^V^^dx^ + — + fdcj) 2 , 
A M dx M = a^V^fidx^ 1 + Adcp, (3.18) 

where the perturbation is specified by the functions {h( T y, h ( T ^, a^r)} of r. As for 
gauge fixing, see Appendix IA.21 The ( LT )- and (T)r-components of the Einstein 
equation give 

h (T)r = (3-19) 
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where C is an arbitrary constant for m 2 = 0, or C = 0 for m 2 7 ^ 0. The (T)-component 
of the Maxwell equation and the (T)^-component of the Einstein equation are reduced 
to 

{r 2 f$l v)i y-V2r 4 A'$l v)2 +m 2 $ iv)1 = 0 , 

2 4 

(^ 2 )' + ^^ + ^ = 0, (3.20) 

where 


$(V )1 

$(y )2 


V / 2a(7’), 

h(T)0 


(3.21) 


For m 2 = 0, if the constant C were nonzero then by going to the gauge h(r) r — 0 
by the gauge transformation (IA.25|) . the metric component Hilt) given by (IA.22|) 
would diverge on the brane. Thus, the regularity of the induced metric requires that 
C = 0 for m 2 = 0. (For m 2 7 ^ 0, it has already been shown that (7 = 0.) 

The regularity of the held strength for the pull-back of Am<Ix m on the brane re¬ 
quires that <E>(vqi is regular at r = r±. With the above equations (13.2011 . the regularity 
of the tetrad components of the six-dimensional Weyl tensor on the brane requires 
that *f)(y )2 / f 1 < h , (y )2 and y/J^yy should be finite at r = r±. These regularity condi¬ 
tions are enough to make equations hermite. Therefore, m 2 is real. 

As in the case of tensor perturbations, we need to rewrite the above equations 
firm into a form which includes only a and which makes it possible for us to take 
the a —> 1 limit without any divergence. Actually, while the variables dWu and 
< h(v )2 were useful to show the reality of m 2 , for numerical calculation there is a more 
convenient choice of variables. The (T)(^-component of the metric perturbation h(r)tp 
in the coordinate (f, ip) is related to h^T)^ as 

h iT )<p = A 6 ^/r + r_/3_h {T)lf . (3.22) 


The coefficient /?_ in the above equation vanishes in the a —> 1 limit. Thus, we rescale 
h(r) <t> as 

hm = (3 ’ 23) 

which approaches to y/A^hyr^ hi the a —> 1 limit. Using this variable, we can rewrite 
the above equations gZDP as 


d 2 


a (T) 


dj , W- 


+ 


f P-? + P+ 


drd(T) + 



/ ((3-r + /3+) 
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2(3- 


(3-r + (3+ 


h(T )< j> ) + 


m 


rj2 1 ! 2(3- 

0= h(T)<b + -- —Orh, 


ftf-r + P+y 
16\/2 


0 -f + /?+ 




i a ( T ) ~ 0 


(/3_f + (3+) V 7i V 5 


rO ( T ) 


-6 


/ 3 - 


/?_f + (3+ 


h(T)<t> + 


?7T 


/ G9_f+ /?+)■ 


<h(T)<j> — 0 , 


(3.24) 


where = Yl f=o 0(1 n as we defined above. The boundary conditions are obtained 
by assuming that $(y)i and <h(y )2 can be expanded in the Taylor series at r = r±: 


dr&(T) + 


8 ^A 6 




2 (3- 


= 0 


r —>±1 


{P-r + (3 +) 4 dj V 7i V 5 V~' ^ (3_f + (3+ 

2 


h 


T)</> 


+ 


m 


{f3- r + (3+) d f f 


a (T) 


= 0 . 


r —»±1 


(3.25) 


3.2.2 Analytic solution for a — 1 


The system of the perturbation equations (EDI becomes simple in the a —► 1 limit, 
and then can be solved analytically. Here we summarize those solutions, which are 
used when we numerically solve the equations (H21 and ■ rm . for general a. Taking 
ol —> 1 and using the equation (Eud, the equations (EOH) become 


d f [(1 — r~)df(i(T)\ + 2A e dfh(r)(p + /J?a(r) 
12 '• - ' /<2 r h, 


- \/ + - _ _ 2 «(7> 


0 , 

0 , 


(3.26) 


where 


h- 


~ 9 

m 

2 A^’ 


(3.27) 


The absence of the zero mode solution can be easily shown using the above differential 
equations and the regularity of the variables at the boundaries EE1 In the following, 
we consider the case of p 2 ^ 0. 

By the change of variables 


Tl = d(T), 

fk 2 = — a/2A 6dfh(j"),p + 2 a,(T)i (3.28) 
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the above equations become 

dr[(l-f 2 )dr^ 1 ]~^2-2^ 1 )+^ 1 = 0 , 

dr[(l-f 2 )dr^ 2 ]+li 2 (^ 2 -^ 1 ) = 0 . 

The original variables are written in terms of T i and To as 

a (T) = Ti = ^2 {®r [(1 - r 2 )dr^ 2 ] + ^2} , 


h 


1 — r 


^2 


<F)V 




f* 2 - 


V 2 A 6 h 2 

This set of equations can be rewritten as 

dr [(1 - f 2 )d f E ± \ + A ±E± = 0, 

where 


E± = dr[(l-f 2 )d^2] +\ t ^2, 

A-t = /i 2 + 1 i y/2/J'? + 1. 

Thus, the general solution is 

^2 = C+P v+ (f) + C_P u _{f) + D + Q v+ (r) + D_Q v _(f), 
where C± and D± are constants and v± is a solution to 

z/ ± (z/ ± + 1) = A±. 


(3.29) 

(3.30) 

(3.31) 

(3.32) 

(3.33) 

(3.34) 


From the regularity of ^(vqi, v7^(v)i> ®(v) 2 // and ^[v) 2 i we can s h° w that 'f'l, 
a /1 — f 2 <9 f T 1 , T 2 and <9 f T 2 should be hnite at the boundaries. At f = 1 , T 2 is regular 
only if D + + D_ = 0 whereas the regularity of is reduced through the differential 
equations to (A+ — /i 2 ) D + + (A_ — /r 2 ) D_ = 0. Since A + 7 ^ A_ for /i 2 7 ^ —1/2, we 
obtain = IA_ = 0. At f = — 1 , the regularity conditions for the variables become 


C .|_ sin u + n + C- sin z/_ 7 r = 0, 

C + A + sin u + tt + C_ A_ sin z/_ 7 r = 0. (3.35) 


Non-trivial solutions can exist if u + e Z or z/_ e Z. In general we can choose non¬ 
negative v±s so that these conditions are explicitly written down as 


v± 


+ 5 ± 4y/2fJL 2 +1. 


= 0 , 1 , 2 , 


(3.36) 


Therefore the KK mass spectrum for vector perturbation is obtained. The case of 
ji 2 7 ^ — 1/2 can be excluded J2BJ- 
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3.2.3 Numerical solution of KK modes 


Here we obtain the first few KK modes of vector type perturbations by numerically 
solving the perturbed Einstein equation (ECU and the junction condition (M - We 
rewrite the system of two second order differential equations (ECT to a system of 
four first order differential equations by defining <9 f a(T) and dfh(T)<f> as well as O(t) 
and h(T)<j> as dependent variables. Here, we take the number of points on a mesh 
M = 101. We solve the problem while changing a from 1 to 0, each time with a 
slightly different value of a. Starging from a = 1, the analytic solutions for a = 1 
presented in the previous subsection are used as a trial solution. 

Figure 01 and 0 show the first four KK mode solutions of a^r) and h for a = 1.0 
and 0.31. The normalization is determined by using the generalized Klein-Gordon 
norm as in tensor perturbations. For vector perturbations, it is defined by 


~. M 4 A</> 

($,V)kg = \ 


d 3 x / drr\ 




QiAviv - + j 

(3.37) 


See the Appendix IA.41 for the derivation. In terms of a(T) and h( T)<j>, this is written as 


($A)KG=(k 0 + k' 0 )6 3 ( k 


k n M$ A0 
’ 2 


dr 


2a[” 1 ) ) aj^ ) 


+ l h 


Ol) / (n-2) 


(3.38) 


where superscript rt*. means that a and are the solutions of vector pertur¬ 
bations with eigenvalue , and we normalized the constant vector as = 1. 

Using the coordinate (f, ip) and the rescaled variable the Klein-Gordon norm 

is further rewritten as 


($,V)kg 


(fco + K) ^ (k 

{ko + k' Q ) 6 3 (k 


k' 


k' 


M 6 4 A<p 
4A 6 

4A fi 


dr 


%(’ ll )/.(" 2 ) I A_h( n l) /J n 2) 
ZU (T) U (T) "T" F U (T)cj> U (T)cl> 


( ni ) T ( rei ) 
U (T) ’ U (T)<t> 


An 2 ) 


a)rP ', h 


f 

(n 2 ) \ 


(3.39) 


We normalize the solution by 


7 GO ui.ni) (na) lfa)\ _ X 
\ U (T ) > n {T)<j) U (T) ’ ,b {T)(j>) °n 1 n 2 - 


(3.40) 


We can easily prove the orthogonality between modes with different m 2 by using the 
equation of motion for vector perturbations. 

Finally, we show the spectrum of m 2 + for the first four KK modes as a function of 
a in figure El We find that m 2 + is non-negative for the entire range of a. Therefore, 
the background spacetime is also dynamically stable in the vector-type sector. 
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Figure 5: The spectrum of m\ for vector perturbations as a function of a. 

The behavior of m 2 + in a 0 limit has the same feature as in tensor perturbations. 
It remains finite in this limit, that is, m 2 + oc a 0 . As was mentioned in subsect ion IT 1 .■'11 
this result is consistent with our previous study |T2j, where we have found that the 
energy scale H* + at which the KK modes begin to modify the effective Friedmann 
equation on the brane at r = r + behaves as H 2 + oc a 0 when a ~ 0. 

3.3 Scalar-type perturbation 

Finally, we show that m 2 > 0 for any non-vanishing scalar perturbations satisfy¬ 
ing relevant boundary conditions. We first derive the perturbed Einstein equations 
and the boundary conditions in subsection ITT II As well as the vector perturba¬ 
tions, there are two physical degrees of freedom. We then rewrite the system of the 
perturbation equations into a form including only the background parameter a. In 
subsection ITT 21 we show that m 2 is real. We summarize the analytic solution for 
a — 1 in subsection ITT Ml Using this result, we numerically solve the perturbation 
equations for a < 1 by relaxation method in subsection ITT41 
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3.3.1 Basic equations 

For scalar perturbations, we can take an analog of the longitudinal gauge: 

(] r 2 

dsl = r 2 (l + 'SfY)r) lJ , v dx* l dx v + [1 + (5>i + «3> 2 )Y] — + Zh^d^Ydx^dct) 

+ [1 - ($! + 3* 2 )Y] /# 2 , 

AMdx M = a r Y dr + (A + a ( f,Y)d(t>, ( 3 - 41 ) 

where perturbations are specified by the functions {T, $ 1 , <f> 2 , h^tp, a r , of r 

and the harmonics Y = exp As for gauge fixing, see Appendix I A. 21 The 

r0-component of the Einstein equation implies that = Cf(r ), where C is an 

arbitrary constant. By using the residual gauge freedom C in (IA.37D . we can set 
C — 0. Thus, 

h(L)<t> = 0 . ( 3 . 42 ) 

The r-component of the Maxwell equation and the (LL)-, ( L)r - and (L)^-components 
of the Einstein equations give 


a r = 0, 

T = $ 2 , 

Aa* = + f'*2, (3.43) 

where a prime denotes derivative with respect to r. The remaining equations are 
reduced to 

0, 

0, (3.44) 

where m 2 = —r]^ u k fl k u . 

With these equations, it is straightforward to show that linear perturbations of R, 
Rm^n'R'mn' 1 K ,ai R-m and R klmn,m ' Rklmn-m 1 are independent linear combinations 
of /$!, (/$i)', $2 and f '&2 — (m 2 /2r 2 )$i and that the matrix made of the coefficients 
remains regular and invertible in the r —> r± limit. Thus, /‘hi, (/ < h 1 ) / , $2 and 
/'$2 — (m 2 /2r 2 ) < f ) i must be regular on the boundaries. 

The correct boundary conditions can be obtained either by using the formalism 
developed in Sendouda et al. 2B] or by setting the coefficients of 1// in (13.441) to zero 


*? + 2( j + ;K 


4A 6 ^ . m 2 

— + $ 2 + -^7^1 = 

/ r 2 f 


4,. m 


2 r 2 / 


($i + 2$ 2 ) = 
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on the boundary, where / vanishes. In the latter method, we obtain 


m 


2/'<h , 1 -4A 6 ($ 1 + <D 2 ) + —$ 


m~ 

2^ 


($! + 2$2 


r—*r± 


r—>r± 


o, 

0. 


(3.45) 


By using the boundary conditions, it is shown that and <h 2 have regular Taylor 
expansion w.r.t. r — r±. In practice, it is useful to note that for given values of $i(r±) 
and $ 2 (r ± ), 


$2 (r±) 
*'i(r±) 

*Z(r±) 






m 

2rf 


(r±) 

f'(r±) ’ 
.2 






$i(r±) 

f'{r±) 


+ 2 <h 2 (r ± ) 


(3.46) 


The above derivation of the boundary condition and the Taylor expansion is simple 
and easy to follow. An alternative and more rigorous derivation is also possible by 
using the formalism developed in Sendouda et al. |2E|- Following the formalism, the 
straightforward calculation gives 


\Ts (\Tff 


’/if) + 

Sr - f ' 2 (v 7 //') V + V r f('/fK*Y - \Sf'(fhrr)' 


r—>7--1- 


r—>r± 


i r^r± 


A'Sr +, 

where h rr = ($! + $ 2 )//, = — ($i + 3<h 2 )/ and & r is defined by 

f Sr + h H = 0. 


By using the regularity of /<hi, (f&i)', <h 2 and /'$' 2 — (m 2 /2r 2 )$! at r = 
boundary conditions are reduced to 


= 0, 

= 0, 
=( 3 . 07 ) 


(3.48) 
r±, the 


/<& i| r=r± = 0, 2/'d> 2 + (/$,)' | r=r± = 0. (3.49) 

This is actually equivalent to (jS33). Indeed, the same Taylor expansion, i.e. (irm 
follows also from firm 
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Heretofore, the system of the above differential equations and the boundary con¬ 
ditions can be rewritten into a form which includes only the parameter a. In terms 
of f, the Ea. (Id.441) becomes 




+ 4 


d f f P- 

J +5 H 


f ' “/3-r + /3 + 

0 - 




m 


/ 


(/?_r + /?+) / 


$1 = 0, 


P-r + P+ 


d f A>, + 


m 


($! + 2$ 2 ) = 0. 
2 (P-r + P + ) 2 f { 


(3.50) 


The function / and thus the Eq. EH3) include only the parameter a. For numerical 
calculations in subsection Id. .'1. 41 we use the first two equations in (Id. 401) as boundary 
conditions for and $ 2 . These boundary conditions are also rewritten in terms of 
f, m and a: 


+ 2$ 2 | f=±1 = 0, <9 f $i 



1 

4 dff 




r==bl 


(3.51) 


Thus, to show the dynamical stability of scalar type perturbations, we calculate the 
spectrum of fh 2 as a function of a by using (ESDI) with (13-511). and show that m 2 is 
non-negative throughout. 

As stated in the beginning of this subsection, we need to follow the four steps (i)- 
(iv) to show the stability. The first step (i) has not yet been considered at all, while 
the set of differential equations ESDI) with ESH) is simple enough and expected to 
be useful for the remaining steps (ii)-(iv). Here we note that the Eci. (Id . 44j) is not a 
manifestly self-adjoint system and, thus, we do not yet know whether the eigenvalue 
m 2 is real or not. In the next subsection ld.d.2l we shall reduce the system of differential 
equations to a manifestly self-adjoint one and show that m 2 is indeed real. After that, 
we shall come back to the equations ESDI) with EHD again and perform numerical 
calculations. 


3.3.2 Reality of m 2 for scalar perturbation 

To show the reality of m 2 for scalar perturbations, we suppose that m 2 ^ 0 and adopt 
the following gauge: 

dv^ 

dsl = r 2 (l + QiY)r] flu dx ti dx u + 2BV( L ' )/1 dx fJ 'dr + (1 + AY)—+ (1 - 3QiY)fd(fi 2 , 

A M dx M = CYdr + (A + Q 2 Y)d4 i, (3.52) 

where Q\. Q i, A, B and C are functions of r. In this gauge the branes are at 
r = r± + 3ff'QiY\ r=r± . Note that the corresponding metric in the five-dimensional 
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Einstein frame after reducing the </ direction is 

dr 2 

ds\( E ) = r 2 f 1/3 r] tlv dx^dx u + 2f 1/3 BV( L ) lx dx lx dr + [1 + (A - Qi)Y]j (3.53) 

Thus, in this gauge a constant-r hypersurface in the five-dimensional Einstein frame 
is flat. Note also that the ( 0 /)-component of the six-dimensional metric and the 
/-component of the U(l) field behave as scalar fields in five-dimension after reducing 
the / direction. Therefore, Q\ and Q 2 represent perturbations of the scalar fields on 
the flat hypersurface in the five-dimensional Einstein frame. 

In the analysis of cosmological perturbations in the four-dimensional Einstein the¬ 
ory, the so called Mnkhanov variables play important roles. The Mukhanov variables 
represent perturbations of the scalar fields on flat hypersurfaces, and thus, are anal¬ 
ogous to our variables Q\ and Q 2 - Therefore, it is expected that the above gauge 
choice simplifies the analysis of perturbations. 

The r-component of the Maxwell equation and the ( LL )- and (L)r-components of 
the Einstein equations give the following algebraic equations. 


C 



A 


0 , 


2 rf , _ 2 rA' _ 12 (A 6 r 2 + 9/)( 2 / + rf) + r 2 /' 2 Qi 

6 f + Gf + rf (6 f + rf’) 2 r 

3(2/ - rf) 4 rA' 

6 f + rf Qf + rf^ 21 


8 A 6 r 3 A' Q 2 
( 6 / + rf 1 ) 2 r 

(3.54) 


where a prime denotes derivative with respect to r and m 2 = —rf w k IJ k ll . The remain¬ 
ing equations are reduced to 

LQ = rn 2 QQ, (3.55) 

where 

«=(£)■ (356) 

and 


L 

ct 

P 

1 


d r ctd r + d r /3 + (3d r + 7 , 
3/ 0 
0 1 


3 4 

-r A 
2 


0 -1 
1 0 

16r 6 A 6 A' 2 
(6 f + rf') : 


24 



f2 


r^(3f 0\ 

/ V o \)' 

9/(2/ + /Q 9(2A 6 + W 2 ) 
rM' 2 32A 6 rM' 2 ’ 

3(2/ - ’7') 

4rA' ■ (3 ' 57) 

With these equations, the linear perturbations of R, R^n'Rmn' ■> R' AI R-,m and 
^■a7n-m'^kl M a f th e positions of the branes are independent linear combinations 
of Q i, Q 2 , Q'i and Q 2 evaluated at r — r± 1 . The matrix made of the coefficients is 
regular and invertible. Thus, Q 1 , Q 2 , Q\ and Q' 2 must be regular at the positions of 
branes. With this regularity condition, it is shown by using the formalism developed 
in ref. [2HJ that the boundary condition at r = r± is 


a = 

b = 


fQl\r=r± ~ Ql\r=r± ~ O' 


(3.58) 


Alternatively, the same bonndary condition can be obtained from the boundary con¬ 
dition and the relation (13.591) below, ft is easy to show that the operator L 

with this boundary condition is hcrmite and that m 2 is real unless Q\ = Q 2 = 0 
everywhere in the interval r_ < r < r + . 

The Mukhanov-type variables {Qi, Q 2 ) can be written in terms of the metric 
variables ($ 1 , $ 2 ) in the analog of the longitudinal gauge. We have shown that m 2 
is real unless the former variables vanish everywhere in the interval r_ < r < r+. On 
the other hand, we have used the latter variables in the numerical calculations. The 
relations between the sets of variables are 


Q 1 
Q2 


$2 + 


+ 


2/$! 

6 f + rf’ 
/riTj 
6 f + rf' 


(3.59) 


where it is understood that a 0 is expressed in terms of $1 and $ 2 - Equivalently, ($ 1 , 
$ 2 ) are written in terms of (Q 1 , Q 2 ) as 

= (6 f + rf)-, 

_ r 

1 Note that the positions of the branes are r = r± + 3f'fQiY\ r - r± and we do not know a priori 
whether fQi vanishes at r = r± or not. Thus, for example, the linear perturbation of R at the 
positions of the branes is SR + 3 f'fQiY evaluated at r = r±, where and SR are the 
background value and the linear perturbation of the Ricci scalar, respectively. Similar statements 
hold also for the linear perturbations of Rm'N'-^mn j R ;M R-,m and f> KLMN '- M R KLMN . M , at the 
positions of the branes. 
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(3.60) 



where it is understood that B is expressed in terms of Q 1 and Q 2 - Therefore, if 
Qi = Q 2 = 0 everywhere in the interval r_ < r < r + then $1 = $2 = 0 everywhere 
in the interval r_ < r < r + . Therefore, we have shown that m 2 is real unless 
$1 = $2 = 0 everywhere in the interval r_ < r < r + . 

3.3.3 Analytic solution for a = 1 

The perturbation equations given in the subsection 13.3.11 can be analytically solved 
for a = 1 as shown in Sendouda et al. eg. Here we summarize the solution obtained 
there. By taking the a —> 1 limit of the equations and (133TD we get 


(l -r 2 ) -4r<9 f «I>i -2($i + <f> 2 ) + fJ?$i = 0, 

(l - f 2 ) <9 2 $ 2 + y ($1 + 24> 2 ) = 0, (3.61) 


and 


$1 + 2 < 3> 2 | f=±1 = 0, (1 - f 2 ) 9 r -Ti| f=±1 = 0, (3.62) 


(3.62) 


where [P = rh 2 / 2A 6 . The solution for the zero mode can be excluded using the first 
condition of (|3.621) and the regularity of the variables at the boundaries (223- In the 
following, we consider the case of /r 2 p 0. 


The differential equations are combined to 


dr[(l - f 2 )d- r F ± \ + X t F ± = 0 , 


(3.63) 


where 


F± = dr[(l - f 2 )a r -T 2 ] + A ± d>2, 


(3.64) 

(3.65) 


\± — / 1 2 + 1 A 3/i 2 + 1. 


If jj 2 p —1/3, the solution is obtained as 


4*2 — C+P v+ + C—P v _ + D+Q u+ + 

$1 = — \{dr[i l ~ h 2 )^f < h2] + 2 fdr§2 + $ 2 }, 


(3.66) 

(3.67) 


where indices v± are determined by 


v±{y± + 1) — A± — jj 2 + 1 i \J 3/i 2 + 1. 


(3.68) 
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The first boundary condition is expanded around r = 1 as 

$1 + 2$ 2 ~ 2/i 2 (.D + + D—)( 1 — r) 1 

-+[(A+ - l)D + + (A_ - 1)D_] 

+0(1-f). (3.69) 

This means D + = D_ = 0, since A + ^ A_ for /i 2 ^ —1/3. Next we expand it around 
f — — 1 , then 

$1 + 2 < 3>2 ~ 47 t _ 1 /i _ 2 (C' + sin u + n + C- sin i/_7r)(l + r ) -1 

+27r~ 1 /W 2 [(A + — 1)^+ sin zi + 7 T + (A_ — 1)C_ sin z/_ 7 r] 

+0(1 + f). (3.70) 

This implies that non-trivial solutions can exist if v + € Z or z/_ G Z. With these 
choices of the parameters, we can confirm that the second boundary condition in (|3.62l) 
is satisfied. In general we can choose non-negative v±s so that these conditions are 
explicitly written down as 

-1 + a/4/i 2 + 5±4v/3/i 2 + 1 

W =- Y -g-= 0,1,2,.... (3.71) 

Therefore the KK mass spectrum for scalar perturbation is obtained. The case of 
/i 2 7 ^ —1/3 can be excluded[28]. 

3.3.4 Numerical solution of KK modes 

Here we obtain the first few KK modes of scalar type perturbations by numerically 
solving the perturbed Einstein equation and the junction condition lism We 

rewrite the system of two second order differential equations (13.5011 to a system of 
four first order differential equations by defining <9 f 4>i and < 9 f < 3>2 as well as 4>i and $2 
as dependent variables. Here, we take M = 51. We solve the problem while changing 
a from 1 to 0, each time with a slightly different value of a. For a — 1, the analytic 
solutions presented in the previous subsection are used as trial solutions. 

Figure E] and [TJshow the first four KK mode solutions of <f>i and $2 for a = 1.0 and 
0.31. The normalization is determined by using the generalized Klein-Gordon norm 
(I A. 6411 given in the Appendix I A. 41 In terms of the coordinate (f, <p), the equation 
dm is written as 

(<M)*g = M 6 4 (2vr) 3 ((u 1 +a;2)5 3 (k 1 -k2)e-^-^^^ 
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dr 


x 


x 


l-i 2 (8/?_(/Lf + P + )d f f + 24/?!/ + A 6 (/?_f + /?+) 2 ) 

+ m*)(/?_f + /?+) 2 $i** + (/ 3 - r + / ?+ ) 4 a _/(^g_ $i + $ 2 a f vi/*) 

(^“ r + ^ ^13p 2 _f + 3A 6 ^~ r + l3+ ^ } 


+(/3_r + /3+) 2 ( 2/?_(/?_ r + /?+)<9 f / + 12/?! / + 3Ag 


(/3 r + /? + ) s 




2^2 


+ (/3 - r +/ + ) 2 ((/3-f + /?+) 2 <9 f / 2 + 10/?_(/?_f + /?+) W + 24/?! P 


4/ 


3A 6 / 


A/3_f + /? + ) 2 


($i + 2 $ 2 )(^i + 2^; 


3(/?f+ /? + ) 3 0 r -/$,tf: 


r=l 


16(8 /?_<9 f / + A 6 (/3_r + /?+)) J f= _ x 
((2tt) 3 ( Wi + cu 2 )5 3 (k 1 - k2 ) e -^i-^)* A ^ r + ? - U u 

ZAq \ 

We normalize the solution by 


= MU 




$ 2 .(3.72) 


$! rai) ,$^ ni) 


®( n2 ) dj( n2 A — A 

V 1 1^2 — °n 1 n 2 t 


(3.73) 


where superscript means that and are the solutions of scalar perturba¬ 
tions with eigenvalue m? nk . 

Finally, we show the spectrum of m 2 + of the first four KK modes as a function of a 
in figure |H1 We find that m 2 + is non-negative for the entire range of a. Therefore, the 
background spacetime is dynamically stable in the scalar-type sector. The behavior 
of m 2 + in a —> 0 limit is also similar to the cases of vector and tensor perturbations. It 
remains finite in this limit, which is consistent with our previous work |Q, where we 
obtained the energy scale at which the correction to the effective Friedmann equation 
on the brane appears. 


4 Summary and Discussion 

In this paper, we have considered the dynamical stability of the six-dimensional brane 
world model with warped flux compactification recently found by the authors. This 
solution captures essential features of the warped flux compactification, including 
warped geometry, compactification, a magnetic flux, and one or two 3-brane(s). For 
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Figure 8: The spectrum of m 2 + for scalar perturbations as a function of a. 

simplicity we have set the four-dimensional cosmological constant to zero, and re¬ 
stricted to linear perturbations with the axisymmetry corresponding to the rotation 
in the two-dimensional bulk. We have expanded perturbations by scalar-, vector- and 
tensor-type harmonics of the four-dimensional Minkowski spacetime and analyzed 
each type separately. 

The perturbations were labeled by its type and values of mass squared m 2 = 
—ri lw k lt k u . To study the stability, we have utilized the fact that for any perturbation 
type, the perturbation equations are reduced to an eigenvalue problem with eigenvalue 
m 2 . We regarded the background spacetime to be dynamically stable if the spectrum 
of m 2 is non-negative. This question was treated in each sector through the following 
steps. 

(i) We have shown the reality of m 2 . 

(ii) We have rewritten the systems of the perturbation equations into a form which 
depends on background parameters through just one parameter a. (0 < a < 1) 

(iii) In the a —> 1 limit we have analytically solved the perturbation equations and 
have shown that the spectrum of m 2 is non-negative. 

(iv) We have numerically evaluated how each eigenvalue m 2 changes as a moves from 
1 to 0, and have shown that the spectrum remains non-negative throughout. 
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For a < 1, we had to numerically solve the systems of the perturbation equations 
many times, each time with a slightly different values of a. In such a situation, 
relaxation method is useful, and thus we have employed this method. When a is 
slightly changed, the previous solution will be a good initial guess, and relaxation 
works well. Since the perturbations can be analytically solved in all the sectors for 
a = 1, we have started the sequence of numerical computations from a = 1 — A a by 
using the analytic solutions as the initial guess, where Aa is a small positive number. 

The mass squared that has physical meaning is = —r^rj^k^k^, which is the 
one observed on the brane at r = r±. We have shown the spectra of m 2 + as a function 
of a in each sector, and found that m 2 + is non-negative for the entire range of a, [0,1]. 
Therefore, the background spacetime we consider is dynamical stable in each sector. 
We have found that there are zero modes only in the tensor sector, corresponding to 
the four-dimensional gravitons. 

Another remarkable feature of the spectra of m 2 + is that they remain finite in a —> 0 
limit, that is, m 2 + oc a 0 and m 2 _ oc a ~ 2 for a —> 0 in all the sectors. This result is 
consistent with our previous study H2. where we analyzed how the Hubble expansion 
rate H± on each brane changes when the brane tension changes. We also considered 
higher-order corrections of the effective Friedmann equation with respect to H±. The 
result is that higher-order corrections appears when H± is larger than a critical value 
H*±. For a —> 0, we found that H*± behave as the equation (13TT1) . Since the higher- 
order corrections are caused by the KK modes, this energy scale corresponds to their 
mass squared. Thus, the behavior of m\ we obtained is consistent with our previous 
result. 

Having established the stability of the exact braneworld solution, there are many 
subjects for future research, including the recovery of 4-dimensional linearized Ein¬ 
stein gravity and corrections to it, the recovery of the 4-dimensional Friedmann equa¬ 
tion, properties of black hole geometries |2HJ ; and so on. 

Note added 

After this paper was submitted for publication, we were notified that the exact solu¬ 
tion considered by the authors had been already found in mum in advance of our 
previous paper [T2j. 
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Appendices 

A.l Harmonics in Minkowski spacetime 

In this appendix we give definitions of scalar, vector and tensor harmonics in an 
n-dimensional Minkowski spacetime. Throughout this appendix, n-dimensional coor¬ 
dinates are aT (p, — 0,1, • • •, n — 1), r)^ is the Minkowski metric, and all indices are 
raised and lowered by the Minkowski metric and its inverse 77 *“'. 

A. 1.1 Scalar harmonics 

The scalar harmonics are given by 

Y = exp (ik p x p ), (A.l) 

by which any function / can be expanded as 

f = I dk cY, (A.2) 

where c is a constant depending on k. Hereafter, k and dk are abbreviations of {AT} 
(ji — 0,1, • • •, n — 1) and n^=o respectively. We omit k in most cases. 


A. 1.2 Vector harmonics 

In general, any vector held v fl can be decomposed as 

W = v (T)n + d„f, (A. 3) 
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where / is a function and V(t)h is a transverse vector field: 


<9 M u (T)/i = 0. (A.4) 

Thus, the vector held v p can be expanded by using the scalar harmonics Y and 
transverse vector harmonics V(t)h as 

= J dk [c(t)V(t)h + C ( L ) 9 M y ] . (A. 5 ) 

Here, cp) and cpq are constants depending on k, and the transverse vector harmonics 
V( 7 > are given by 

V(T)n = u M exp(ifc p a: p ), (A. 6 ) 

where the constant vector satishes the following condition. 

= 0 (A.7) 

for k^kfj, 7 ^ 0 , and 

k^u^ = 0 , 

= 0 (A. 8 ) 

for non-vanishing k fl satisfying k^k^ = 0 , where is an arbitrary constant time¬ 
like vector. For k M = 0, the constant vector tA does not need to satisfy any of the 
above conditions. For the special case = 0, the second condition in (1A.8I) can 
be imposed by redefinition of cpq. Actually this condition is necessary to eliminate 
redundancy. Note that the number of independent vectors satisfying the above con¬ 
dition is n — 1 for k^k^ ^ 0 and n — 2 for k^k^ = 0 and that these numbers are equal 
to the numbers of physical degrees of freedom for massive and massless spin -1 fields 
in n-dimensions, respectively. 

Because of the expansion (lA.Sjl . it is convenient to define longitudinal vector 
harmonics by 

V {L) » = d tJ Y = ikfjY- (A.9) 


A. 1.3 Tensor harmonics 

In general, a symmetric second-rank tensor held t^ v can be decomposed as 

t/u/ T v T d u v p T frilly, (A. 10) 
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where / is a function, v u is a vector field and t(r)nv is a transverse traceless symmetric 
tensor field: 


"( 2 > 

d^t(T)fiv 


0 , 

0 . 


(A. 11) 


Thus, the tensor field t^ can be expanded by using the scalar harmonics Y, the 
vector harmonics V(t) and V(l), and transverse traceless tensor harmonics T(t) as 


t 




j dk \c(T)T(T)iiv + C^LT){d^V{T)v + d^V^T)^) 
+C(LL){d^y(L)v + + C(Y)Ylj^ J,] . 


(A-12) 


Here, C(r), C(lt ), C(ll), and C(y) are constants depending on k , and the transverse 
traceless tensor harmonics Tm are given by 

T(t)^ = exp(ik p x p ), (A.13) 


where the constant symmetric second-rank tensor s^ satishes the following condition. 


for k^kfj, 7 ^ 0 , and 



0 , 

0 


k^s^ 

■V 

T ^ 1 S 
' ° \±v 


o, 

0 , 

0 


(A. 14) 


(A.15) 


for non-vanishing k^ satisfying k^k^ = 0 , where is an arbitrary constant timelike 
vector. For k M = 0, the constant tensor s^ does not need to satisfy any of the 
above conditions. For the special case k^k^ = 0, the last condition in (IA.15I) can be 
imposed by redefinition of c^lt), c (ll) and C(y)- Actually this condition is necessary 
to eliminate redundancy. Note that the number of independent symmetric second- 
rank tensors satisfying the above conditions is (n + l)(n — 2)/2 for k^k^ ^ 0 and 
n(n — 3)/2 for k^k^ = 0 and that these numbers are equal to numbers of physical 
degrees of freedom for massive and massless spin -2 fields in n-dimensions, respectively. 

Because of the expansion (1A. 121) . it is convenient to define tensor harmonics T^lt), 
and T(y) by 

T^lt)^ = d^V{r) u + 
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(A.16) 


i('U/ j./Cj, “l - 

T(LL)fiv = dfiY(L)i/ + dvV^L)/! ~ ~V^9 P V^L) p 

= ^- 2 k p k v + ^-k p k p q, 1 ^j Y, 

T(y)(u/ = q fu/Y. 

A.2 Gauge transformation of the perturbations 

In this appendix we give gauge transformations of the perturbations of the metric 
and the 17(1) gauge field. The coordinate gauge transformation is of the form 

x M ^x M + f M . (A. 17) 

There is another kind of gauge transformation. The perturbations of the field strength 
SFmn are not changed under gauge transformation of the gauge held, 

5A m — 5 - 5A m + <9mC- (A. 18) 

The gauge parameters f M and ( can be expanded by the scalar and the vector har¬ 
monics as 


II 

l-Uy 1 

= {f(T)V(T) P + £(l)V(l) m ) dx p + f r Ydr + f^Ydf, 

(A. 19 ) 

c = a", 


(A. 20 ) 

where fr, ^ and C 

are supposed to depend only on r. Under the above gauge 

transformation, the perturbation variables transform as 



h{T) = h(T), 

(A. 21 ) 


h(LT) = h(LT) - €(T), 

(A. 22 ) 


h(LL) = h(LL) - £(L), 

(A. 23 ) 


h(Y) = h (Y) - 2 rffr + ^k p k p f ( L), 

(A. 24 ) 


2 

h(T)r = h (T )r + ~£(T). ~ C(T)> 

(A. 25 ) 


2 , 

h(L)r = h(L)r + ~£(L) ~ £(L) — 

(A. 26 ) 


k{L)<j> = h(L)<t> f<f>, 

(A. 27 ) 


f' 

h — h - - —f 

lL rr IL rr Sr *5 

(A. 28 ) 
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h — 

i i r (j) 

a (T) = 
a {L) = 


^r(f> 


£ 4 , + y£</» 


h<W - }fir, 

a (T), 

r A t 

a (L) + C — y? 0 i 


a r 


a r ""t“ &rQ j* 



^(f) ^<j) f A £ 7 * * 



(A.29) 

(A.30) 
(A.31) 

(A.32) 

(A.33) 
(A.34) 


Finally we summarize the gauge conditions employed in this paper. Of course, 
tensor type perturbations are gauge invariant from the beginning. For vector pertur¬ 
bations, we set h(LT) = 0 by taking £,(t) = h(LT) ■ The master equations for scalar 
perturbations in the main part of section 13.31 are derived in the analog of the Lon¬ 
gitudinal gauge where h(LL) = h(L)r — h r <j, = 0. This can be obtained by taking a 
gauge 


£(L) = h(LL), 

(A.35) 

2 

£,r ~ h(L} r + ~h(LL) ^(LL) > 

(A.36) 

^-Cf(r) + f{r) J dr' , 

(A.37) 


where C is an arbitrary constant representing the residual gauge freedom. We use 
this residual gauge freedom to set C = 0 in subsection 13.3.11 In subsection 13.3.21 
we show the reality of m 2 for scalar perturbations by adopting another gauge. This 
gauge corresponds to h(LL) = 0 , 3 fh(y) + r 2 h ( j >( j > = 0 and h r( j, = 0 , which can be set by 

£(L) = h(LL), 

f _ 3 fh( Y ) + r 2 h H + 3 k^k^ LL) 
rf ( 6 / + rf) 

it = Cf(r ) + f{r) J dr 'yyy, 

where C is again an arbitrary constant representing the residual gauge freedom. For 
the U(l) gauge held, we set a^) = 0 by ( = both in the main part of section l3~3l 
and in subsection 13.3.21 


(A.38) 
(A.39) 

(A.40) 
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A.3 Relaxation method 


In the following we explain relaxation method in detail [M]. First of all, we rewrite 
a system of second-order differential equations to a system of first-order differential 
equations of the form 



where yi denotes one of N dependent functions. For example, y\ = h and 1/2 = d?h for 
tensor perturbations. When we numerically solve the equations, one of the dependent 
functions has to be fixed at some f in an arbitrary manner. (When we plot the 
solution, we use another normalization.) Thus, these dependent functions have to 
satisfy N + 1 boundary conditions instead of just N. The problem is overdetermined 
and in general there is no solution for arbitrary values of m 2 . For certain special values 
of m 2 , the Eq. (|A.41 1) does have a solution. Such rh 2 is the eigenvalue. It is convenient 
to reduce this problem to the standard case by introducing a new dependent variable 


Vn+i = m 2 


(A.42) 


and another differential equation 


dyN+i = q 
dr 


(A.43) 


Next, we replace the differential equations with finite-difference equations. We 
first define a mesh by a set of k — 1, 2, • • •, M points at which we supply values for 
the independent variable f k - I 11 particular, f\ (= —1) is the initial boundary, and 
Vm (=1) is the final boundary. We use the notation y*, to refer to the entire set 
of dependent variables yi, 2 / 2 , • • ■, Vn+i at point fk ■ At an arbitrary point k in the 
middle of the mesh, we approximate the set of five differential equations by algebraic 
relations of the form 

0 = E fc = y k - y fc _i - (f k - f fc _i) g fc (f fc , r k - 1 , y k, Yk~ 1 ) , k = 2, 3, • ■ ■, M. (A.44) 

The finite-difference equations labeled by E fc provide a total of (N + 1) (M — 1) equa¬ 
tions for the (N + 1 ) M unknowns. The remaining N + 1 equations come from the 
boundary conditions. At the first boundary we have 


0 = Ei = B (r 1, yi) 


(A.45) 


while at the second boundary 


0 — Em+i = C (tmi Ym) ■ 


(A.46) 
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The vectors B and C have only ri 2 and n\ — N + 1 — ri 2 nonzero components respec¬ 
tively, corresponding to number of the boundary conditions. 

The solution of the Eq. (I A. 441) . (IA.45I) . and (IA.46j) consists of a set of variables y*, 
at the M points f k . The algorithm we now describe requires an initial guess for y*.. We 
then determine increments Ay*, such that y^ + Ay^ is an improved approximation to 
the solution. Equations for the increments are developed by expanding the Ea. (IA.44j) 
in first-order Taylor series with respect to Ay^. At an interior point, k — 2, 3, • • •, M 
this gives 


0 = E k (y fc + Ay k , y fc _i + Ay*,^) ~ E k (y fe , y fc _i) 

N+l N +1 


E 


A„. , A n . 

~^yn,k —1 “ 1 “ / J ^Vn^ki 


/ r\ ‘— m £jn,K— 1 ' / r\ ‘—*£//<.,/t - } (A.47) 

n=l ’ „=1 ^4 

where y n , k is the value of y n at the point r k . This provides (N + 1) (M — 1) equations 
for Ay k . Similarly, the algebraic relations at the boundaries can be expanded in a 
first-order Taylor series for increments that improve the solution. Since Ei depends 
only on y 1; we find at the first boundary: 

n+i (9E 

0 = E, + V (A.48) 

A d V' 


n.l 


At the second boundary, 


0 = E 


M+l 


N+l 

\ ^ 9E m +i a 
/ , —-A y n , M . 


(A.49) 


—7 dy n ,M 

We again note that the Ea. (IA.48j) and (IA.49I) have only n 2 and n k = N + 1 — n 2 
nonzero components. We thus have a set of (N + 1) M linear equations to be solved 
for the corrections Ay*,, iterating until the corrections are sufficiently small. 


A.4 The generalized Klein-Gordon norm 

To normalize mode functions we use the generalized Klcin-Gordon norm, whose def¬ 
inition can be easily read off from the kinetic term in the effective action for the 
corresponding physical degree of freedom. See, for example, Appendix A of ref. [32] 
for the definition and the motivation of the generalized Klcin-Gordon norm. 


A.4.1 Tensor perturbation 

For tensor perturbation, 

(Jr *^ 

ds 2 & = r 2 [q, JV + h^] dx^dx" + — + /d</> 2 , 
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(A.50) 


A M dx M = Adqb, 

where h pv is a symmetric, transverse and traceless four-dimensional tensor depending 
on (x^, r), the bulk action is expanded up to the second order in perturbation as 

A = ^yf d 6 xV^ (R - 2A 6 - ^F mn F mn ^J 

= ~ I dAx J drr 2 ‘r] pp '‘ri aa ' [lf^df.hpvdvhp'v’ + f r 2 d r h pa d r h p ' a '](A.5l) 

Here, we have not written down the boundary term since it does not change the 
definition of the generalized Klein-Gordon norm. It is easy to check that the correct 
equation of motion is derived from this action. 

From this form of the action we can read off the generalized Klein-Gordon norm 
as 

($, F) kg = -iMiM. I rf 3 x f drr 2 r ji i i 1 ' r j vu ' ■ (A.52) 


A.4.2 Vector perturbation 


For vector perturbation, after fixing the gauge freedom (h^r) = 0) and using the 
corresponding constraint equation (the ( LT ) component of the Einstein equation), 
the metric and the U( 1) held in the linearized level are written as 


dsl = r 2 r]p U dx p dx l/ + 2 
A M dx M = andx p + Ad(f>, 


, dr + h^dcj) 

L r f 


dv^ 

dx p + — + fd(j .) 2 , 


(A.53) 


where a p and h,^ are transverse four-dimensional vectors depending on (x p , r ) and 
h r/1 is a transverse four-dimensional vector depending only on x p . The bulk action is 
expanded up to the second order in perturbation as 


r _ M 6 

h ~ ~T 


M|A0 


J d 6 xVA (r - 2A 6 - i F mn F) 

d 4 x [ drL , 


(A.54) 


where 


L = -r] pa 


2rf v dp i a p d ll a (J + —7) ,LV d ll hpd l ,h rJ + 2fr 2 d r a p d r a a — Ar 2 A'hpd r a a + r 2 d r h p d r h a + 6 h p h, 


f 


(A.55) 
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We have not written down the boundary term since it does not change the definition 
of the generalized Klcin-Gordon norm. It is easy to check that the correct equations 
of motion are derived from this action. 

From this form of the action we can read off the generalized Klcin-Gordon norm 
as 


($,V) KG = -iMeM f rf 3 x f drrj 


JIV 


*iA**u, - i,) + - ($ 2 An, - 


f 


2 


)2v 


(A.56) 

where the solutions <f> and T are specified by (*hi M , $ 2 //) and (\h lM , respectively, 

as 


4> : 

T : 


&I1 


K 

K 

V2 7 


hp = r 2 <f> 2m , 
h/x = r 2 T 2/1 . 


(A.57) 


A.4.3 Scalar perturbation 


For scalar perturbation, after fixing the gauge freedom (h^L) = h{L)r = h(L)<j> = 
h r< p = O(l) = 0) and using the corresponding constraint equations (the (LL), ( L)r , 
( L)(j) and ref) components of the Einstein equation and the (L) component of the 
Maxwell equation), the metric and the U(l) field in the linearized level are written 
as 


ds 2 6 = r 2 ( 1 + < h 2 )h ' flv dx fl dx v + 1 + ($1 + $ 2 ) 


dr 2 

T 


+ 


1 - ($1 + 3<f>2 


A M dx M 


- < A+ * 


d 5 (/r 2 * 1 )' + /'i 2 


d(p, 


fd4 > 2 , 


(A.58) 


where $1 and $2 are functions of (x M , r). The bulk action is expanded up to the 
second order in perturbation as 


where 


h 


^ J d 6 xVA (r - 2A 6 - ^F MN F MN ^j 



L = -rTd^ndu Q + Q r LQ 


(A.59) 

(A.60) 
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Here, L and ft are defined in ()3.57j) and 


Q 


(Qim 

V &[*] 



(A.61) 


where 


Qim 

Q2$] 


$2 + 


2 /d>! 


1 

~A' 


6 f + rf n 

1 


2 r 2 


(/r 2 <f>i)' + /'$ 2 


+ 


frA'&i 
6 / + rf r 


(A.62) 


We have not written down the boundary term since it does not change the definition 
of the generalized Klcin-Gordon norm. It is easy to check that the correct equations 
of motion are derived from this action. 

From this form of the action we can read off the generalized Klcin-Gordon norm 
as 


rr+ 


/ crx / drr z 


x 


(<3i[*]$<3;[*] - e;[*]ft(ji[*]) + - - Owe 


where the solutions $ and are specified by pairs of five-dimensional functions ($i, 
$ 2 ) and ('hi, d' 2 ), respectively. For = <Fj(r)e* J? ' il/fc i xI/ and 4/ = dh(r)e M? ' i! ' fc by 
using the equations of motion, integrating by part and using f(r±) = 0 , we obtain 


= M 6 4 A0(27r) 3 (a; 1 +a;2)5 3 (k 1 -k2)e-^-^ t 


rr+ 


dr 


r _ 8(2r/' + 6/ + A 6 r 2 ) 


x 


- (m 2 + ml)r 2 ^i^l + 2 r 4 / , ('k 2 < 9 r < f>i + < £ , 2 < 9 r'I / i) 


+r 2 (13/ + 3A 6 r 2 ) < F 1 'h^ + 4r 2 (2rf + 12/ + 3A 6 r 2 )4> 2 ^ 

2 

+ j(r 2 f 2 + 10 rff + 24 f 2 + 3A 6 r 2 /)(<F 1 + 2$ 2 )(^ + 2^*) 


3r 3 / / $i'hJ 


Ll6(2/' + A 6 r)J r=r _( ’ 


r=r+ 


(A.64) 


where Ui = k® and m 2 = It is shown that ($, ^)kg — 0 for m 2 7 ^ m 2 and 

that ($, ^)kg is time independent. 
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